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Abstract 

A Hartree-Bose mean-field approximation for the IBM-3 is presented. A 
Hartree-Bose transformation from spherical to deformed bosons with charge- 
dependent parameters is proposed which allows bosonic pair correlations and 
includes higher angular momentum bosons. The formalism contains previ- 
ously proposed IBM-2 and IBM-3 intrinsic states as particular limits. 
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With the advent of radioactive nuclear beam (RNB) facilities unexplored regions of 
the nuclear chart will become available for spectroscopic studies. New aspects of nuclear 
dynamics and novel types of collectivity and nuclear topologies are expected. Nuclei with 
roughly equal numbers of protons and neutrons (Z ~ N) and with masses in between 
40 Ca and 100 Sn are of particular interest because the build-up of nuclear collectivity in this 
mass region occurs in the presence of pairing correlations between alike nucleons (proton- 
proton and neutron-neutron) as well as between neutrons and protons Q . This offers the 
possibility of experimentally accessing nuclei that exhibit a superconducting phase arising 
from proton-neutron Cooper pairs || . Although recent breakthroughs |3],f|] have made shell- 
model calculations possible for this mass region, they are still of a daunting complexity and 
alternative approximation schemes are required that yield a better intuitive (e.g. geometric) 
insight. 

One of the possible alternatives is the Interacting Boson Model (IBM) ||. It has been 
shown || that nuclei with protons and neutrons filling the same valence shell require an 
extended boson model, IBM-3. In IBM-3 three types of bosons are included: proton-proton 
(tt), neutron-neutron (is), and proton-neutron (5). The tt, is, and S bosons are the three 
members of a T — 1 triplet, and their inclusion is necessary to obtain an isospin-invariant 
formulation of the IBM. Over the last decade the validity of the IBM-3 has been tested and 
its relationship with the shell model worked out |7j]- [PH . 

The mean-field formalism has been an important tool to acquire a geometric under- 
standing of the IBM ground state and of the vibrations around the deformed equilibrium 
shape fll2|- |14}| . Moreover, a treatment based on mean-field techniques generally leads to 



a considerable reduction in the complexity of the calculation, allowing the introduction of 
additional degrees of freedom if needed. Studies in the intrinsic framework are thus useful 
to assess the importance of higher angular momentum bosons with e.g. £ = 3~, 4 + , . . . or to 
investigate the role of extra degrees of freedom not included in IBM such as two-quasiparticle 
excitations, etc. 

An intrinsic-state formalism for the IBM-3 was recently presented by Ginocchio and 



Leviatan (GL) [|nj. In that work charge-independent deformation parameters are imposed 
in the Hartree-Bose transformation from spherical to axially deformed bosons and the trial 
wavefunction is taken to have good isospin. Closer inspection reveals that this trial wave- 
function has the additional isospin SU(3) symmetry which in IBM-3 is equivalent to orbital 
U(6) symmetry. (Isospin SU(3) symmetry is to IBM-3 what F-spin symmetry ]TB is to 



IBM-2.) Moreover, it is well known that isospin symmetry itself is increasingly broken in 
Z ~ N nuclei as the nuclear mass increases [ 17| . There is also tentative evidence for rigid 



triaxial shapes in the region of interest and its proper description would require the inclusion 
of three body forces or higher angular momentum bosons. We therefore present in this pa- 
per a generalization of the treatment of GL in which none of the above symmetries (isospin 
SU(3) and SU(2)) is imposed on the trial wavefunction and which includes bosons of angular 
momenta higher than £ = 2. For practical applications we restrict ourselves here to £ = 0, 2. 

We start with the usual spherical boson creation and annihilation operators 7j mr , jemr 
where £ is the angular momentum, m is its third component, and r is the isospin projection. 
Each boson carries isospin T = 1. We also define 7^ mT = (— l)^~ m 7f_ mT - In terms of these 
boson operators, a system of N bosons interacting through a general number-conserving 
two-body Hamiltonian can be written in multipolar form as 



tf = J> r7 ] T -7^ + E E <T 2 r,rA T2 -fL TA , (1) 



J . z,, _l r L ^ L - ^ L 
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where the symbol ■ denotes scalar product in orbital space. In isospin space the only re- 



striction is r% + r 2 = r 3 + r 4 (i.e., a charge- conserving Hamiltonian is assumed) and T^ T2 are 
multipole operators with total angular momentum L, 

^M,rir 2 = E yUxliwSCIixTx X 7^ 2 t 2 )m) (2) 
his 

where the coupling is only done in angular momentum. The Hamiltonian ([I]) can be used 
for IBM-3, for IBM-2, or even for a general isospin non-conserving Hamiltonian with three 
kinds of bosons. 

Deformed bosons are defined in terms of spherical ones by means of a unitary Hartree- 
Bose transformation 



Em P 

and their hermitian conjugates. The deformation parameters rf^ in these equations verify 
the orthonormalization conditions 

%m Vim = V> 1^ Vim Vi'm' = <W<W- (4) 

Im P 

Note the explicit dependence on the isospin component r of the transformation 77, allowing 
different structures for the different condensed bosons ir, is, and 5. The index p labels 
different possible deformed bosons. We choose p = for the fundamental deformed bosons 
and p — 1,2,... for the different excited bosons. For instance, in an SU(3) scheme different 
values of p = 0, 1, 2, 3 label the ground, /3, 7, and scissors bands, respectively. Since in this 
work we only treat the ground-state condensed boson, the Hartree superscript p is always 
zero here and it will be omitted in the following. The formalism for the excited states will 
be presented elsewhere. 

Following Ref. [15|j, the trial wavefunction for the ground state of an even-even system 



with a proton excess is of the form (the trial wavefunction for an even-even system with a 
neutron excess is obtained by interchanging the role of protons and neutrons) 

10(a)) = A^( a )T\ Np - Nn \0), (5) 

where the operator creates a correlated bosonic pair in isospin space 

At(a) = T\VU + aVlTl (6) 

In Eq. ([5]) N p (N n ) is the number of proton (neutron) pairs in the valence space. The trial 
wavefunction (0) contains the isospin-conserving formalism of GL and the IBM-2 as natural 
limits. Two different values of a are connected with these limits. For a = — |, A' (a) 
corresponds, in the particular case of r-independent deformation parameters, to an isoscalar 
bosonic pair. Its total isospin is T — N p — N n and the results of GL are reproduced. Any 
other value of a breaks isospin symmetry. In particular, a = eliminates the mixing of 5 
bosons in the ground state and yields an IBM-2 intrinsic state. It should be emphasized 



that when the deformation parameters rf^ in Eq. (3) depend on the isospin component r, 
the set of operators with r = —1, 0, 1 do not form an isospin triplet, and consequently 
the At in Eq. (6) may contain mixtures of T=0,l,2 isospin components. 

At this point we would like to remark that the trial wavefunction (|5]) is not the most 
general U(18) intrinsic state, involving a combination of all r = — 1, 0, 1 condense bosons. 
The U(18) intrinsic state is written as 

l^(i 8 ) = (rt) Np+Nn \o), (7) 

where 

i^E^wyLr- (8) 

imr 

In this state orbital angular momemtum, isospin, and charge are broken. The state given in 
Eq. (HP improves over this state by including charge conserving pair correlations. Thus, it 
is expected to lead to a deeper energy minimum. Numerical results illustrating this point 
will be presented later on. 

The variational parameters of the trial wavefunction are the matrix elements rf irn of the 
Hartree-Bose transformation, associated with the orbital and isospin degrees of freedom, 
and the parameter a, which determines the amount of mixing of 5 bosons in the ground 
state. 

The ground-state energy is obtained by taking the expectation value of the Hamiltonian 
(H|) in the state (§): 

E(r], a) = J2 £rfi(a, r) + £ V T c uT2 ^ T J 2 (a, n^n), (9) 

T T1T2T3T4 

where 

e r = Ew&L ( 10 ) 

im 

\rc T/ *Tx *T2 T3 T4 / -i -i \ 

n,T"2,T3,T4 — 2^ * £imiTi,l2m2T2,i3'm3T3,£4m4T4 % 1 mi%2m2^3m3^4m4' l-L-U 
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(0(g)| rtiy 
(0(g) | 0(g)) 

and 

v (0(^)1^^^^ |0(g)) 

/a g, riT 2 r 3 r4 = 1 2 . 13 

(0(g) | 0(g)) 

The coefficients Eg T include the single particle energies eg T in Eq. (]T[) plus contributions from 
the two body term in the same equation. The coefficients Ve imiTlt £ 2m2T2t i 3m3T3t i 4m4T4 are the 
symmtrized interaction matrix elements between normalized two-boson states following Ref. 

Ve imi T 1 ,e 2 m2T2,e 3 m 3 T3,e 4 m 4 T 4 = \ (hmiTi, h^^l H l^m^Ts, t^m^r^) 

(14) 



The dependence of the energy on the variational parameters 77's is contained in the one- 
body e (10) and the two-body V c (|i"T|) terms, while the dependence on a comes through the 



isospin matrix elements f\ (|12D and / 2 (13). The latter matrix elements are straightforward 
to calculate by a binomial expansion of the ground-state trial wavefunction (|5|). 

The Hartree-Bose equations for the orbital variational parameters rj are obtained by 
minimizing the energy (|]) constrained by the norm of the transformation. Assuming a 
charge- conserving Hamiltonian ([I]) the following Hartree-Bose equations result: 

^Lxmx,t 2 m 2 r fl 2 m 2 = ^rVtimi' (^) 

l 2 m 2 

where the Hartree-Bose matrix h T is 

(16) 

lim 3 l 4 ra 4 T2T-j,T 4 <2 m 2 

The term i]^ 2m2 in the denominator is a consequence of a mathematical trick for obtaining 
a set of three coupled Hartree-Bose equations fljjl) . These depend on the isospin matrices 



fi and f2- For each value of a the matrices fi and f'2 are calculated and the Hartree- 
Bose equations (p~5HT6|) are solved self-consistently. The procedure is iterated until one finds 



the absolute minimum of the energy (|9|). Once the problem is solved self-consistently, the 
diagonalization of (|15D provides the deformation parameters r}\ m for the ground state. 

To test the present formalism and to compare with the one by GL, we used a simple 



Hamiltonian recently proposed by Ginocchio [nj 



H=-k £ P T :P T , (17) 

T=0,l,2 

where 

P T =( s i~d+(-l) T S~s) L=2 > T . (18) 

In these equations the symbol : denotes a scalar product in orbital and isospin spaces and 
Jimr = (—^y~ m+1 ~ T Je~m-T- The Hamiltonian ([17]) is clearly isospin invariant and provides 
a first simple test to the present formalism 

Figure 1 shows, for a system with 5 proton pairs and 3 neutron pairs, the ground-state 
energy for the Hamiltonian QT7D as a function of a. The dashed line is calculated with 
r-independent deformation parameters; the GL minimum energy is reproduced for a = — |. 
The full line is calculated with the present formalism. The latter calculation always gives a 
lower energy and, in particular, the minimum is not obtained for a = — | but for a ~ —0.32. 
In addition, the corresponding deformation parameters are r dependent. The energy gained 
by breaking isospin invariance in our trial wavefunction is relatively small. In this respect 
it may be advantageus to use the GL intrisic state for isospin conserving hamiltonians. 
Though, a better approximation would be obtained by performing variation after isospin 
projection over our trial wavef unctions. 

We note that for a system with equal number of protons and neutrons the present 
formalism recovers exactly the GL results; differences occur for Z 7^ N. This can be seen 
in Fig. 2 where the deformation parameters (3 T are plotted versus the difference N p — N n 
(starting with 4 proton pairs and 4 neutron pairs). The deformation parameters (3 T are 



obtained from (3 T = J , ^ T , 2 — T (see Eq. (1) of Ref. [|T^]). For N p = N n the deformation 



l* T l 

parameters are independent of r, but not any longer as N p — N n increases. The proton and 
neutron deformations remain very close; the 5 deformation /3g, however, quickly becomes 
very large in comparison. This is because N$ decreases as N p — N n increases. This effect 
can be seen in Fig. 3 where the mean values of the boson numbers, iV T , are plotted. The 
same behaviour has been obtained recently with large scale shell model calculations (see 
Ref. ||20|). In all our calculations we found that the Ginocchio Hamiltonian ( |TTD leads to a 
7-independent energy surface. 

It is worth noting that the present formalism allows one to reproduce the well-known 
case of triaxiality in IBM-2. To show this we use the IBM-2 Hamiltonian 

H = -(Q^ + Q' U )-(Q 7T + Q' U ) , (19) 

where • denotes scalar product in angular momentum, Q v is the SU(3) generator, Q = 
s^d + d)s — ^{d) x d) L=2 , for proton bosons and Q' u is the SU(3) generator, Q' = s^d + 
efts + 2 x ^) L=2 ' f° r neutron bosons. The minimization procedure now gives a = 0, 
which corresponds to the IBM-2 limit. In addition, the minimum deformation parameters 
correspond to a prolate proton condensate, axially symmetric about the intrinsic z axis, and 
to an oblate neutron condensate axially symmetric about the intrinsic y axis, giving rise 
to an overall triaxial shape. It should be pointed out that in this case there is no triaxial 
minimum for aligned proton-neutron shapes with equal deformations. Here the overall 
shape is triaxial but the underlying separate proton and neutron condensates correspond to 
different (prolate-oblate) axial shapes. 

Finally, we present a calculation in which isospin is explicitely broken by the hamiltonian: 

# = £ e tT n iT - \ N[Q° : Q° + ^Q 1 : Q 1 ] (20) 

It b 6 

where Q T = + d^s — ^(cft x d)] L=2,T and N[...] stands for normal ordering product. In 
Figs. 4 and 5 we present the results of a calculation with e S7T = e su = 0, e d7T = e du = 1.5, 
e sS = 2.3 \N P -N n \ and e d5 = 1.5 + 2.3 \N P + N n \ (all e's in MeV). In Fig. 4 the deformation 
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parameters, /3's, are shown as a funtion of N p — N n (the calculation starts with N p = N n = 4 
and then N p is increased). Fig. 5 shows the corresponding ground state energies for the 
intrinsic states of GL, and thos defined in eq. 5 and 7. It is interesting to note that 
the GL intrinsic state and our pair correlated intrinsic state produce the same results for 
N p — N n = 0, being better than the U(18) intrinsic state. For N p — N n > both isospin 
nonconserving intrinsic states give better results than GL. It is also interesting to see that 
our pair correlated ansatz is superior to the the U(18) for moderate values of N p — N n . No 
triaxial deformation is found in these calculations. 

In summary, we have extended the intrinsic-state formalism of Ginocchio and Leviatan 
|15|| for IBM-3 in three different ways. First, the Hartree-Bose transformation is chosen to 
depend on the isospin component r. Second, variable isospin bosonic pair correlations are 
introduced through the parameter a. Finally, higher-order bosons, other than the usual s 
and d bosons, are included in the Hartree-Bose transformation. This formalism contains the 
IBM-2 and GL intrinsic states as particular limits. Substantial differences in the deformation 
parameters are obtained when N p ^ N n . We have presented results for isospin conserving 
and nonconserving hamiltonians with s and d bosons. Substantial differences in the defor- 
mation parameters j3 T are obtained for N p > N n . In most of the cases studied, our pair 
correlated intrinsic ground states are lower in energy than the GL ground states, although 
for the isospin conserving hamiltonian the energy gain is small. We therefore conclude that 
the new intrisic state is useful for treating isospin breaking hamiltonians. 

This work has been supported in part by the Spanish DGICYT under contracts No. 
PB95/0123 and PB95-0533, a DGICYT-IN2P3 agreement and by the European Commission 
under contract CIl*-CT94-0072. 
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FIGURES 

FIG. 1. Calculated ground-state intrinsic energy as a function of a for a system with 5 proton 



and 3 neutron pairs interacting through the Ginocchio Hamiltonian (17) with k = 1 MeV 



FIG. 2. Deformation parameters j3 T for a system with N n = 4 neutron pairs as a function of the 
difference N p — N n between the numbers of proton and neutron pairs. The Ginocchio Hamiltonian 
(|i~7|) with k = 1 MeV is used. 



FIG. 3. Mean values of the boson numbers, N T , for a system with N n = 4 neutron pairs as a 
function of the difference N p — N n between the numbers of proton and neutron pairs. The Ginocchio 
Hamiltonian ( |Tt[) with k = 1 MeV is used. 



FIG. 4. Same as Fig. 2 but for the non-conserving isospin Hamiltonian (|20| ) with the parameters 
given in the text. 

FIG. 5. Calculated ground-state intrinsic energy, as a function of the difference N p —N n between 
the numbers of proton and neutron pairs, for the non-conserving isospin Hamiltonian ( |20| ) with 
the parameters given in the text. 
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